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An optimal path to transition in a duct 

By Damien Biau and Alessandro Bottaro 
Die AT, University of Genova, 1 via Montallegro, Genova, Italy 

This paper is concerned with the transition of the laminar flow in a duct of square 
cross-section. Like in the similar case of the pipe flow, the motion is linearly stable 
for all Reynolds numbers, rendering this flow a suitable c andidate for a stu dy of the 
'bypass' path to turbulence. It has already been shown (JBiau et al.ll2008[) that the 



classical linear optimal perturbation problem, yielding optimal disturbances in the 

f— ^ ' form of longitudinal vortices, fails to provide an 'optimal' path to turbulence, i.e. 

fN| , optimal perturbations do not elicit a significant nonlinear response from the flow. 

Previous simulations have also indicated that a pair of travelling waves generates 

immediately, by nonlinear quadratic interactions, an unstable mean flow distortion, 

responsible for rapid breakdown. By the use of functions quantifying the sensitivity 

of the motion to deviations in the base flow, the 'optimal' travelling wave associated 

to its specific defect is found by a variational approach. This optimal solution is 

M \ then integrated in time and shown to display a qualitative similarity to the so- 

^^' called 'minimal defect', for the same parameters. Finally, numerical simulations 

of a 'edge state' are conducted, to identify an unstable solution which mediates 



^ ■ laminar-turbulent transition and relate it to results of the optimisation procedure. 

M-H , Keywords: duct flow; optimal perturbations; minimal defects; transition to 

(^ ' turbulence 

o 

^' 1. Introduction 

, MH' Transition to turbulence in ducts i s still a challenging issue despite the 125 years 



since the seminal paper by Osborne lRevnol ds (1883) which led to the definition of a 

dimensionless number, that came some time afterwards to bear his name, capable of 

K* \ broadly discriminating cases of 'streamlined' flow from cases where 'sinuous motion' 

^^ ' prevailed. Osborne Reynolds recognized also that there was no unique value of 

^— N \ this dimensionless parameter, representing the ratio of convective to viscous forces, 

f^ ' separating the two classes of motion, and that the end state was influenced by 



the background perturbations present. The problem of laminar-turbulent transition 
^^ ■ in shear flows was posed, to fascinate and attract the attention of thousands of 

f^ \ researchers in the years to come. Incidentally, the story of how the paper was 

reviewed by two referees of t he caliber of Lord Raylei gh and Sir George Stokes 

makes for instructive readi ng (I Jackson &: LaundenbOOTI) . 



In a subsequent paper, iRevnoldsl (|1895[ ) goes beyond transition and presents 



for the first time, the decomposition of the turbulent field into a mean and a 
jrt ' fluctuating part, arriving at the equations now known as the 'Reynolds-averaged 

Navier-Stokes equations', with unknown turbulent stress terms in the mean flow 
equations. Reynolds had the intuition that a balance is necessary to maintain tur- 
bulence, namely production, i.e. the transfer of kinetic energy from the mean flow 
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to the fluctuations, must equilibrate dissipation of the fluctuations, for turbulence 
to subsist. 

While some classical flows, such as that between differentially heated parallel 
plates in the gravity field or that between concentric differentially rotating cylinders, 
exhibit a smooth progression to increasingly complicated patterns via a sequence 
of bifurcations, whose initial points are well predicted by (modal) linear stabil- 
ity the ory, the laminar square duct and the pip e flows are linearly asymptotically 
stable ( Gilllll965l iTatsumi fc YoshimuralllQQOl ). On the other hand, experiments 
clearly show that these flows undergo transition abruptly at moderate values of 
the Reynolds number. The transition to turbulence originates from a subcritical 
instability which requires perturbations of finite amplitude to bring the flow out 
of the basin of attraction of the laminar state. Current understanding ascribes the 
failure of classical stability theory to its focus uniquely on the least stable eigen- 
mode. When a small disturbance composed by a weighted combination of linear 
eigenfunctions is considered, because of the non-normality of the linearized sta- 
bility operator, there is the potential for very large transient amplification of the 
disturbance energy, even in nominally stable flow cond i tions . This property has 
been beautifully described in a paper bv lTrefethen et al.l (|l993[ ). where the authors 
also show that non-normal operators render the modal instability of pseudo-modes 
(i.e. modes of a perturbed operator) possible. The definition of pseudospectrum ha s 
later been extended to the case of base flow uncertainties bv iBottaro et al.l (20031). 

A recent new direction in rationalising the abrupt transition of many wall- 
bounded shear flows focuses on identifying alternative solutions (beyond the lami- 
nar state) to the governing Navier-Stokes equations. In the last two decades such 
coherent solutions, generally unstable and in the fo rm of steady states or travel- 
ling waves, have been found, first for channel flows (jNagatalllQQOl . IClever fc Busse 
1992 . Cherhabili fc EhrensteinI 1997 . Wal cfFc 1997|) and more recently in a pipe 
(JFaisst fc Eckhardtll2003i IWedin fc Kerswell ,200411 . The solutions computed have 
the same basic structure, which appears to be in relation to the self-sustaining pro- 
cess (SSP) observed in the near wall region of turbulent shear flows. The three main 
ingredients of the SSP are: streamwise rolls of amplitude 1/Re which induce streaks 
(an order one spanwise modulation of the ideal laminar flow) by lift-up effect and, 
to close the loop, a travelling wave (an unstable eigenmode of the streak) whose 
quadratic interactions feed onto the rolls. 

These alternative solutions to the Navier-Stokes equations, also labelled 'exact 
coherent structures', have been compu ted by continuation m ethods. First the rolls 
are generated by a thermal body force ( Clever fc Busselll992), by a cent rifugal force 
(|Nagatalll990l) . or by an unphysical ad hoc body force (|Waleffelll997|) . The solu- 
tion is then obtained by progressively decreasing the forcing imposed to create the 
rolls. Flow states similar to those computed by continuation approaches have been 
observed in numeri cal simulations of turb ulent flows. By reducing the size of the 
computational box. iHamilton et al.l ( 19951 ) have found the minimal flow unit, con- 
taining just a single pair of streaks, capable of maintaining some turbulent activity. 
Periodic solutions have also been used to describe the bursting event, i.e. the break- 
up a nd re-creation of coherent structures such as streaks, in the turbulent buffer 
laver. lKawahara fc Kidal (|200lh have found two periodic solutions in plane Couette 
flow which represent quiescent and turbulent phases of the flow. The phenomenon 
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of bursting has thus been interpreted in phase space as the wandering of the flow 
trajectory back and forth between such solutions. 

Tracking the bifurcation of SSP states by continuation with the Reynolds num- 
ber has led to two kinds of solutions: the so-called upper and lower branch states. 
From a dynamical system point of view, the lower branch solutions sit on a sep- 
aratrix, the phase-space boundary between the basin of attraction of the laminar 
flow and that of the turbulence. The lower branch flow is also called an 'edge state'. 
An open question concerns the number of different states embedded in the separa- 
trix. These exact coherent solutions, with their stable and unstable manifolds, are 
typically low-dimensional saddle points which collectively produce a chaotic repel- 
lor. This leads to the perspective of representing the chaotic dynamics of turbu- 
lence through a projection onto a (hopefully not too large) set of 'exact' coherent 
solutions. The idea is that these structures, connected through their stable and 
unstable manifolds, are capable to support chaotic trajectories in phase space. A 
chaotic behaviour with a finite number of degrees of freedoni could be explained 
by the existence of a strange attractor (JRuelle fc Takenjl 19711) . The appealing as- 



pect is that turbulence could be described by a finite number of coherent solutions, 
which constitute the skeleton around which the chaotic dynamics is organised. How- 
ever, the question of whether a global attractor of the Navier-Stokes equations in 
three-dimensions exists is an issue that is not yet settled. Beyond its fundamental 
character, this question has important practical implications, for knowing that a 
finite dimensional attractor exists would guarantee that the long-time behavior of 
the Navier-Stokes solutions can be approximated by numerical means (i.e. using a 
finite number of degrees of freedom) . 

The present paper deals with the flow through a duct of square cross-section with 
constant pressure gradient dPo/dx in the streamwise direction. While sharing with 
the cylindrical pipe flow many characteristics (both flows are linearly stable, for ex- 
ample, and undergo spontaneous transition at comparable values of the Reynolds 
number), the square cross-section, by its geometric features, has the capacity of 
strongly constraining secondary flows, both instantaneous and time-averaged. An 
example of this is provided in figure [T] (left frame) where turbulent mean veloc- 
ity components arising from a direct numerical simulation at Re — Urh/v — 150 
in a duct of (dimensional) length equal to Gtt/i are displayed, and the presence 
of symmetric secondary vortices of Prandtl's second kind is shown. The carte- 
sian coordinates employed in the following are a;, y and z to define, respectively, 
the streamwise, vertical and spanwise directions. The scales used to normalise the 
Navier-Stokes equations are the friction velocity Ur = y^{—h/Ap)dPo/dx, the side 
length of the square /i, the density p and the kinematic viscosity v. 

Figure [T] (right), obtained from a series of direct numerical simulations at pro- 
gressively varying values of the pressure gradient, shows that in the turbulent regime 
the Reynolds number Rcb = Ubh/v, Ub bulk velocity, grows almost linearly with Re 
(noted Rb j in this figure f or imm ediate distinction from Reb), in agreement with 
results by lUhlmann et al.l ( 20071 ). The smallest Re for which turbulence exists is 



found to be 130, corresponding to i^Cb » 1730. This value is sligh tly smaller than 
that i dentified bv lUhhnann et al.l (|2007l ) (Re = 160, Reb = 2154). IUhlmann et al " 



(|2007t ) have interpreted Re = 160 as the value defining a 'minimal channel' ca- 
pable of some form of turbulent activity, with two streaks, each flanked by a pair 
of quasi-streamwise vortices, active over two walls which face one another. Such a 
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Figure 1. Left: turbulent mean crossflow vortices and streamwise fiow contours. Spacing 
between adjacent isolines is equal to 4 in units of Ut. Right: variation of Rei, with Re^ in 
the laminar and turbulent regimes. 



state with two pairs of rolls adjacent to opposing walls, with the active role of the 
walls alternating in time, exists near m arginal condition s, and it has been inter- 
preted as the signature of the edge state ( Biau et al.ll2008[ ). Interestingly, the effect 
of the long time averaging of the two possible sets of marginal states (with the two 
pairs of rolls lying on the horizontal or the vertical walls) is that of producing a 
eight- vortex pattern such as that displayed in figure [1] (left). 

The very first step of the bypass transition process in boundary layers, channels 
and ducts, has often been described as an algebraic instability of the base flow, 
leading to the formation of streamwise-elongated streaks. Subsequent steps in the 
process are the linear wavy instability of the streaks and the nonlinear feedback 
onto the rolls. Despite the fact that streaks, rolls and waves are all present in the 
exact coherent solutions briefly cited above, there is as yet no description of how 
such states might arise from the stable laminar base flow. A path to transition 
relying on the amplification of so-called optimal disturbances and their subsequent 
nonlinear development has proven to be sub-optimal, when viewe d from the point o f 
view of the minimal initial energy needed to provoke transition ( Biau et al.ll2008l) . 

A different optimization strategy is thus called for, to identify the optimal path 
from the laminar to the turbulent state, relying on recent developments on the SSP 
and on physical understanding of the breakdown process. 



2. The optimisation approach 

Evidence collected indicates that key to transition is the establishment of a distorted 
mean flow profile, capable of supporting strong instabilities. This distortion is dif- 
ferent from that caused by the rolls and the streaks issued fr om classical optimal 
perturbation theory. The recent theory of 'min imal defects' ( Bottaro et al.l 12003 . 
Biau fc Bottaroll2003. iBen-Dov fc Cohenll2007Qfl l. yielding the base flow deforma- 
tion of fixed amplitude capable of maximal destabilisation, is also not sufficient to 
explain available observations, since it provides no link between the initial state 
and the deformed flow. 

Thus, we have set up a new optimisation strategy, aimed at finding the travelling 
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wave of minimal amplitude capable of yielding a self-sustained state, based on the 
following initial steps: 

• linear algebraic growth of the travelling wave, 

• generation of weak streamwise vortices by quadratic interactions, 

• production of a strong streak by lift-up, 

• regeneration of the travelling wave, closing the loop for a sustained state. 

Velocity and pressure are decomposed into a base, laminar state {Uo,Po), and 
a time-dependent part. The time-dependent component is itself decomposed into a 
slowly time- varying mean flow defect (Q) plus a streamwise travelling wave (q) of 
order 0{e) <^ 1. The vortex {V, W) is of order 0{e^), while the streak U resulting 
from lift-up is O(e^Re). Using the decomposition: 



Uo{y,z) ' 




- U{y,z,t) - 




' u{x,y,z,t) 
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Viy,z,t) 
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v{x,y,z,t) 





W{y,z,t) 


w{x, y, z, t) 


Po{^) . 




. P{y,z,t) 




. p{x,y,z,t) 



the equations governing the mean, streamwise-averaged flow read: 



Vy+W, = 0, 



Ut + ViUa + U)y + W{Uo + U), 

Vt 
Wt 



Re 



[Uyy + Uzz) - {Vu\y + WU | ^ ) , 



-Py + J^iVyy + V,,) - {vv\y + WVI) , 

-Pz + -^{Wyy +Wzz) - {Vw\y +Ww\z) , 



associated to homogeneous boundary conditions U = V^W^Qon the walls. 
Note that in the equations above overbars over the products of fluctuating velocity 
components denote averaging along the streamwise distance. The Navier-Stokes 
equations, linearized around the streaky flow (C/q + C/), are: 



Ux+Vy +Wz = 0, 

Ut + {Uo + U)ux+v{Uo + U)y + w{UQ + U)z = -Px 

Vt + {Un + U)Vx = -Py 

Wt + (Uq + U)wx = -pz - 



Re 
1 

Re 
1 

Re 



{Ux 
{Vx: 
[Wx 



U. 



■Uz 



yy 

yy + Vzz), 



w. 



yy 



■Wzz), 



together with u = v = w = Oon the walls and streamwise periodicity. 

The energies of mean flow defect (E) and fluctuations (e), to be used later, are 
defined as: 



Pi*) = ^ / (U^ + V^ + W^) dy dz, 

2 Jyz 
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e{t) — — {u + V + w'') dx dy dz. 



2 



xyz 



A gain is defined as G = e(i)/e(0), and the goal of tlie work is to find tfie 
initial wave at i = capable of maximising G. In the first step of the iteration 
procedure, U, V and W are zero; they are created from the second iteration on, 
because of the presence of the Reynolds stress terms in the mean flow equations. 
As an alternative it might have been interesting to test the approach with a different 
objective functional, such as a gain based more directly on the Reynolds stresses 
or the ratio between production and dissipation terms. 

The optimisation procedure is standard and, instead of carrying out a con- 
strained optimisation, we introduce a Lagrangian functional C to be maximised 
without constraints. This augmented functional is: 

£ = G-< Q\ L{Q, q)>-< q\ 1{Q, q) >, 

where < . > denotes the usual inner product (i.e. volume integration of the product 
of the two quantities over the domain) and Q\ q'^ are Lagrange multipliers. On an 
extrcmum point, the following necessary conditions must be satisfied: 

d£ dC 

— ,5Q = and ^5q = 0, 

oQ oq 

leading to the adjoint problems in symbolic form: 

< L^Q\ SQ>=- < q\ l{SQ)q >=< Gq, SQ >, 

and 

< /^q^ 5q>= - < Q\ L{5q)Q > = < gq, 5q >, 

(superscript f designs adjoint quantities) which must be solved together with the 
direct equations for the fiuctuations and the mean fiow; SQ and Sq are infinitesimal 
variations of the mean flow and of the fluctuations, Gq and gq are sensitivity 
functions which arise naturally after integration by parts when isolating 6Q and 

Sq. 

In extended form, the adjoint equations for the mean flow are: 

vj + w! = 0, 

-Ul = l-{uly + ul) + Gu, 

^w} + uHuo + u).. = -pi + j-iwl + wl), 
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with U^ — V^ — W^ = on the waUs. The adjoint equations for the fluctuations 
are: 

= 0, 





ul + vl + wl 




-ul - (Uo + U)ul 


^l 


- (Uo + U)vl + u^Uo + U)y 


4- 


-iUo + U)wl+u\Uo + U), 


ut 


= tjt — lyt — on the w 
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pl 
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Re^'^- 
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yy 


+ 4J 


+ 9v, 


pI 
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i?e^"^- 
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""1% 
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w\^ 


) + 9 



with uT — ^t — ^T — on the waUs and periodic boundary conditions along 
X. Note the minus signs in front of the time derivatives of the 'adjoint momentum 
equations', indicating that the only possible direction of stable evolution is negative 
time. All the adjoint equations are linear but are coupled to the direct state, and 
this represents a numerical challenge since the direct fields must be stored at all 
time steps. 

The sensitivity terms are: 

Gij = -u^Ux + Vx +Wx)x + iu^v)y + {u^w)z, 
9u = -U^v - Ulw, 
9v = -Ulu~2Vy^v + w{V,} + W^), 
9^ = -Ulu-2Wlw + v{V^^+Wl). 
Gjj corresponds to the s ensitivity to va r iation s in the base flow; it is a generalisation 



of the function found bv lBottaro et al.l (|2003l ). The terms gu, 9v and g^ correspond 



to wave sensitivity functions, and are necessary for the optimisation of the feedback 
loop. 

A sequence of direct and adjoint calculations is carried out to optimise G by 
employin g spectral collocat ion codes which are slight modifications of those em- 



empioym g spectral collocat ion codes wnicn are slignt modmcations ot tnose em- 
ployed by iBiau et al.l ( 20081 ). The result of the integration of the direct system up 



iot = T (where T is the target time of the optimisation) provides q{T) (and Q{T)); 
the adjoint system is integrated backward in time and is initialised by (l\T) — q{T) 
and Q^(T) ~ 0. Once the adjoint fields are computed, the direct calculations are 
re-initialised with g(0) = g^(0) and Q(0) — 0, so that the cycle can restart. At 
convergence we obtain the optimal initial fluctuating field q{0) whose energy is 
maximised at t — T, and the associated base fiow distortion Q{t) produced by the 
action of the fluctuations onto the mean field. 

The disturbances are expressed using a single-mode Fourier decomposition in 
the streamwise direction, i.e.: 

q{x, y, z, t) = qiy, z, i)e+'"^ + q* (y, z, i)e-'"^ 

where superscript * denotes complex conjugate. The present representation is ac- 
ceptable if the first stages of the transition process are focussed upon. 

3. Optimisation results 

In the following, the Reynolds number is fixed at the value Re = 150. Figure [H 
shows results of the optimisation procedure for the case a = 1, eg = 3 x 10"'^, and 
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for a final target time T = 1. Initially, the transient behaviour of the fluctuation 
energy follows the linear curve (shown with a dotted line) obtained in the absence 
of a base flow defect. Within the optimal procedure developed here the energy E 
of the mean flow distortion increases, as a result of Reynolds stress feedback. For 
t larger than about 0.8 the deformed mean flow can sustain a rapid amplification 
of the fluctuations, which brings the disturbance energy e to a level unattainable 
by the transient growth process alone. Thus it appears that transient amplification 
and defects coiicur in defining the initial stages of transition, as postulated by 
Biau k Bott"arol ()2004l ). 




Figure 2. Behaviour in time of the energy e of the optimal initial fluctuations and of the 
energy E of the optimal defect. The dotted line is the result of a simulation of the linear 
disturbance equations in the absence of mean flow distortion, starting from the optimal 
initial disturbance field. 



Figure [3] reports results from a detailed parametric study of the problem, under 
the conditions indicated in the caption of the figure. It should be preliminarily be 
observed that gains of order 800 can be a chieved by the t ransient amplification 
mechanism in the linear case when a = ( Biau et al.ll2008[) . However, for a = 
no mean flow distortion is produced in a nonlinear setting, and at large times the 
disturbances simply decay. For a = 1 the results are more interesting. 

The two frames labelled as (a) show that for T sufficiently long both e and E 
can increase very much. Within the time interval considered, the energy gain of the 
fluctuations is optimal for T « 1.8, i.e. over a physical viscous time scale of order 
H/ut, which is the characteristic time of evolution of the streaks. At the imposed 
value of eo, the defect plays a role only for T larger than about 0.5; the result 
of optimisations runs carried out for smaller target times coincide with classical 
optimal perturbation results (dashed curve). 

The effect of varying eo is displayed in frames (b) ; as a result of nonlinearities 
the gain can increase very much, causing subcritical transition when the initial 
disturbance energy level is sufficient. 

The effect of the computational box size is represented in frames (c). There are 
two noteworthy results here. The first is that there appears to be a clear cut-off 
length L\. w 250, below which the cycle is not self-sustained. This is the 'minimal 
box dimension' for transition to appear, and the result is consistent with direct sini- 
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Figure 3. Parametric study, (a.) a — 1; eo = 3 x 10"'^. (b) a = 1; T = 1. (c) 
T = 1; eo = 10~^. The dashed lines represent the hnear case, in the absence of mean 

flow defect. 



ulations bv lUhlmann et al.l (|2007f) . The second interesting point is that, for T = 1, 
there is an optimal box length, equal to Lj, ~ 350 (the corresponding wavenumber is 
a = 2'KRe/ L\. « 2.7), where maximum amplification of both fluctuations and mean 
fiow defect is achieved. These results are mitigated by the observation that a single 
wavenumber disturbance has been considered in the model, as already pointed out 
earlier. In reality, one should expect streamwise-Fourier modes interactions when 
long computational boxes are numerically computed. When a single mode with a 
very small is considered (i.e. L|, larger than 2500) the results obtained here are 
superposed to the (classical) optimal case; the amplification process coincides with 
the phenomenon of lift-up of low speed streaks from the wall, with a base flow which 
is not distorted from the ideal con ditions. This demonstr ates that, although clas- 
sical optimal perturbation theory ( Butler fc Farrelllll992r) yields initial fiow states 
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capable of large growth for small a, these are inefficient at triggering the sequence 
of processes which eventually causes breakdown of the flow. The reason is that no 
defect can be generated by the Reynolds stress terms, as shown in frame (c) (right) 
of figure |3l The optimisation strategy outlined in this paper, while recovering re- 
sults of the classical theory for a small, highlights the importance of closing the 
feedback loop. 

4. The minimal defect 

The previous section demonstrates the importance of mean flow defects for the 
sustainement of fluctuations. It is thus interesting here to introduce a link be- 
tween optimal perturbations and minimal defects, i.e. those base flo w distortions of 
minim al energy capable to induce efficiently subcritical instability (|Bottaro et al. 



20031 ). For transitional flows (i.e. low amplitude mean flow distortions), the concept 



of minimal defect provides a new outlook on the possibility of rapid disturbance 
growth. For turbulent flows (distortions of flnite amplitude), these defects display 
characteristics of possible relevance to the dynamics of coherent eddies. 

The following is a brief summary of the methodology to compute minimal de- 
fects. Let us consider the linear stability equations in symbolic form as an eigenvalue 
problem: 

L{Uo)q = -iujq. 

A steady base flow defect 5U induces variations on the perturbations and on the 
corresponding eigenvalues, i.e. 



t/o + 6U 



u + Suj 
q + Sq 



The linear stability problem then becomes: 

L{Uo + 5U){q + 5q) = -i{uj + Suj){q + Sq), 
the first variation of which is 

L{Uo + 6U)q + L{Uo)Sq = —iuiSq — iSujq. 
In order to isolate the variation of u, we use the adjoint variable g^ defined by: 
< gt^ {L + iuj)q >^< q, {L^ ~ iuj*)q'' >= 0, 



with an appropriate inner product ( Bottaro et al.l 120031) . We can thus obtain the 
sensitivity function Gjj'. 

Suj =i < q'f, L{Uo + SU)q > 
= < Gu, 6U > 

with the normalization < g^, g >= 1. Note that here the function Gij differs from 
that obtained earlier. The minimal defect is found by maximising the growth rate 
u!i — Imag{uj), under the constraint that the energy of the defect J {U — Uq)'^ dy dz 
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is fixed at the value E/^ . This latter constraint is a simplified state equation for the 
base flow deviation. The Lagrangian functional takes the form: 

C = uji-x( f {U -UqY dydz - Ea) , 

and, by imposing that (dC/dU) SU — 0, the optimisation loop is: 




ImagiG^jY dydz/AEA 
= Imag{G1j)/2\ 



with the exponent n denoting the iteration number. An illustrative result is pre- 
sented in figure S] (left frame). The qualitative similarity between the results in 
the left frame and that on the central frame (corresponding to the optimisation 
of section 2) is clear, with a low speed streak at the center of wall fianked by two 
high speed streaks, although the minimal defect displays much stronger gradients 
of the velocity and presents a much larger (local) instability growth rate. It has to 
be stressed the fact that the minimal defect is a local feature which does not satisfy 
the momentum conservation equation. When it is imposed as initial condition in 
a direct simulation, the sharp gradients diffuse rapidly leading to a solution which 
resembles more to that displayed in the two right frames. The figure in the right 
frame of figure |3] is very similar to that arising from the optimisation procedure, 
aside from the smaller peak absolute values of U attained. In the direct numerical 
simulation on the right frame of figure 2] nonlinearities contribute to a redistribu- 
tion of the energy among many flow harmonics, whereas in the central frame the 
only wavenumbers which can accommodate the available energy are a — Q and 
a = 1. Despite this, the agreement is sufficiently good to provide confidence in the 
ability of our simple model to provide correct flow features in the initial phases of 
transition. It is also interesting to observe that t he flow patterns disp l ayed i n figure 
m resemble the asymmetric solution computed bv lPringle fc Kerswelll ( 20071 ) for the 
case of the pipe. 

An open question concerns the mech anism which associates the position of each 
streamwise vortex to its sign of rotation. lUhlmann et al.l ( 2007 ) propose a heuristic 



explanation based on a simple kinematic analysis of the interaction of streamwise 
vorticity with a corner. The interaction of a streamwise vortex, close to one of the 
corners, with the impermeable wall can be modelled using three image vortices and 
potential theory. Thus the vortex would migrate because of the induced velocity 
field towards a stable position. This argument is consistent with observations. Con- 
sidering the results in figure 0] another hypothesis can be formulated: a low speed 
streak sitting on a bisector is the 'most active' fiow feature, i.e. it is that particular 
flow structure capable of better sustaining the wall cycle. 

If we express the wavy perturbations as q(a;, y, z, t) = A{t)q(y, z)exp{iax — iujrt) , 
the variation in time of the amplitude, when a mean flow is mildly distorted, obeys 
the equation: 

— — - = uji + Suji = uji+ < Imag{Gu), SU > . 
A dt 

The product < Imag{Gu), SU > can be interpreted as a Landau coefficient which 
determines the nonlinear character of the bifurcation; if this coefficient is positive 
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Figure 4. Left frame: the minimal defect is plotted with isolines of the streamwise dis- 
turbance velocity, for Ea = 0.6598 and a = 1. Center frame: snapshot of the streak U 
obtained by the optimisation procedure of section 2. The figure corresponds to t = 0.6 in 
figure[21 Right frame: instantaneous result (i = 0.6) of a direct simulation of the transition 
process (corresponding to the continuous black line of figures [5] and [T]) initialised with the 
optimal solution of figure [51 In all three frames the spacing between the U isolines is equal 
to one (i.e. Ut in dimensional terms), and the dotted lines denote negative values. 



then the nonlinearity is destabilising and the bifurcation is subcritical. It is inter- 
esting here to note that 5U is always proportional to A^Re whether we consider 
the lift-up regenerating loop (through the formation of streamwise vortices) or the 
direct generation of streaks by the Reynolds stresses in the streamwise momentum 
equation (and the latter assertion is clear after realising that the characteristic time 
scale of evolution of the defect is the viscous scale). It is thus feasible to assume that 
a second self-sustaining mechanism exists, with the streaks directly re-generated by 
the fluctuations. 



5. Nonlinear simulations of the edge state 

Numerical simulations of non-trivial states are described in this section; the calcu- 
lations at Re =150 are initialised at t = with the laminar flow plus the optimal 
wavy perturbation of figure 2 (with a = 1 and eo = 3 x 10"'^) multiplied by an 
appropriate scaling factor /3. The time evolution of the skin friction is depicted in 
figure [5] and the blue trajectory which relaminarises after a few units of time has 
/3 = 0.92. The edge of chaos is a stable manifold residing in phase space somewhere 
between the lamina r and the turbulent flow s. To find this invariant object the pro- 
cedure described by ISchneider et al.l (|2007l) has been adopted, successively refining 
(every viscous time unit) the initial guesses on either sides of the edge surface in 
phase space, through a continuous update of the relaxation parameter (3. Four inter- 
mediate solutions (a couple initiated at time t = and a couple initiated at i = 15, 
on both sides of the edge) are also shown with colours in figure [5l highlighting the 
exponential separation with t within each pair, with the trace of one state swinging 
up to turbulence and the other relaminarising. 

The edge state lives in between and is represented by the dashed line. Snap- 
shots of the secondary flow patterns along the dashed curve are provided in figure 
[SI It is interesting to observe that the basic building block of the edge state is 
represented by two pairs of vortices. The larger outer pair is weaker, it sits above 
a smaller near-wall pair, which presents an upwash region fluctuating around a bi- 
sector. The patterns displayed in the figure are similar to those shown in figure 4; 
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Figure 5. Skin friction versus time. Dotted lines indicate the laminar and the (mean) 
turbulent values of the skin friction. The dashed curve corresponds to the edge state. Blue 
trajectories relaminarise, the black and the red trajectories lead to turbulence (before 
eventual relaminarisation at large times). 



here, as time varies, different walls become 'active'. The time-averaged edge state 
in a pipe of circular cross-section consists of one strong pair of vortices near the 
surface with a much weaker p air above it ([Schneider et al.l 120071 lEckhardt et al. 
2007llPringle fc Kerswel]ll2n07h . 



t = 0.6 



t = 33 



t = 40 



t = 54 







Figure 6. Instantaneous secondary flows of the solution on the edge. 

In turbulent shear flows, the spatio-temporal chaos is extensive, i.e. the number 
of degrees of freedom scales with the system's volume. Thus it is impossible to 
repres ent a turbule n t flow as a dynamical system in a low dimensional phase dia- 
gram. iToh fc Itanol (J2003ll used the total energy input versus the dissipation rate 
as a two-dimensional picture of the dynamics. Here, in figure [71 we use a reduced 
phase diagram representation spanned by the energy of the streamwise-averaged 
flow, Eu = \ J„ ^(U + Uo)^ + V'^ + W'^ dydz, and by the bulk Reynolds number (see 



also Biau et al.[ 



The black curve, i.e. the trajectory that goes to turbulence shortly after t — 
0, shows that chaos is attained very rapidly (on a convective time scale) once 
the trace starts diverging from the edge. After a few orbits around an unstable 
(turbulent) node the flow is ejected along the unstable-laminar manifold. During 
the relaminarisation process, the patterns of motion remain self-similar while slowly 



Article submitted to Royal Society 



14 



D. Biau and A. Bottaro 



350 




1800 2000 2200 2400 2600 2800 3000 3200 
Re, 




2900 



Figure 7. Eu 'energy' versus Rcb- The laminar fixed point with Reb = 3163 and 
Eu = 306.45 is denoted by the square symbol. The figure on the right frame provides 
better details of the fiow trajectory on the edge (dashed curve). 



decaying. The x— structure disappears first and the flow tends to be formed by 
quasi-straight rolls and streaks, which then decay on a viscous time scale towards 
the laminar fixed point. The dashed curve in figure [7] shows the trajectory which 
oscillates on the edge surface: this state is a relative attractor and all trajectories 
which are initially on the edge surface remain confined to it. Conversely, if the initial 
condition is close to the edge surface, but not exactly onto it (the case of the red 
and blue curves), the trajectory will arrive close to the saddle, before being ejected 
away. 

It is significant that the optimal initial condition obtained from the procedure 
outlined in section 2 yields a solution which resembles that found on the edge and 
follows closely, for some time, the dashed curve in figure [5l before departing away 
from it. This fact signals that the simplified optimisation procedure that we have 
devised produces an appreciably good estimate of the 'true optimal' (as also attested 
by the fact that the relaxation parameter /3 needed to maintain the solution on the 
edge is close to one), where by 'true optimal' we define that fluctuating state of 
minimal initial energy cq capable to provoke transition. The initial disturbance of 
minimal energy causing breakdown arrives arbitrarily close to the hyperbolic point, 
and, of course, the closer it gets, the more slowly transition is triggered. 

6. Conclusions 

Transition in many wall-bounded shear flows arises from a complicated interaction 
of rolls, streaks and waves; 125 years ago Osborne Reynolds was able to catch a 
glimpse of this interaction. Efforts since, to unravel the secrets of transition to 
turbulence, have focused mostly on linear stability theory (modal and nonmodal) 
and, more recently, on 'exact coherent states', the self-sustaining process and, to a 
lesser extent, on 'minimal defects'. We have tried to put together these individual 
bricks to provide a coherent picture of the process. 

The new optimisation approach described here provides for the first time a com- 
plete model of a transitional path including, in particular, the feedback of the fiuc- 
tuations onto the mean fiow. We have shown that the initial stages of transition rely 
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on a combination of algebraic and exponential amplification of disturbances, with 
the latter closely associated to the creation of a mean flow defect. If disturbances 
of long streamwise wavelength are initially excited by the receptivity conditions, 
mean flow distortions cannot be created, except in the obvious case of very noisy 
environment. 

Albeit simplified, our optimal model of transition has produced a solution which 
sits initially on a trajectory directed towards the hyperbolic point on the edge 
surface, i.e. the separatrix between the laminar basin of attraction and the chaotic 
dynamics. Direct simulations have confirmed the suitability of the model proposed. 
The basic flow structure of the edge state in the cross-section of the square duct is 
formed by two overlapping pairs of vortices, with the near-wall pair more intense. 
This pattern resembles that obtained by minimal defect th e ory a n d is also similar 
to the asymmetric edge sta t e coin pute d by Schneider et al.' ('20071) . Eckhardt et al 



(|2007l ). |Pringle fc Kerswelll ()2007l ) and lPuguet et all (12008) for the pipe flow. 



Work in progress focuses on the effects of increasing the Reynolds number, to 
try and compute the unsteady coherent solutions encountered during turbulent 
intermittency or in puffs. From a practical point of view, it should be possible to 
use the optimisation technique introduced here also to determine semi-empirical 
transition crite ria based on a si mple balance between production and dissipation, 
as indicated bv iRevnolda (|l895[ ). 



The support provided by the EU under Marie Curie grant EST FLUBIO 20228-2006 is 
gratefully acknowledged. 
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